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We provide a short summary of three-flavour neutrino oscillation parameters as determined from 
a global fit to 2008 data, and we comment on possible hints in favour of a non-zero value of the 
mixing angle found in arXiv:0806.2649. We do confirm a hint from solar + KamLAND data 
at about 1.5(7, which can be understood from the recent SNO CC/NC measurment. However, 
we show that a claimed hint from atmospheric data is much less robust. It depends on details 
of event rate calculations and treatment of theoretical uncertainties. We could identify two data 
points showing an "excess" (at the 1 o level) in the SK-I multi-Ge V e-like data, which seem to 
be the origin of the slight preference for 613 > 0. We point out that once SK-I and SK-II data are 
combined this "excess" disappears, and irrespective of the details of the analysis, no "hint" from 
atmospheric data is obtained for the SK-I and SK-II combined data set. As a result the global fit 
of all data leads to a best fit value of 613 consistent with zero within less than Icr. 
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1. Introduction 

Thanks to the synergy amongst a variety of experiments involving solar and atmospheric neu- 
trinos, as well as man-made neutrinos at nuclear power plants and accelerators neutrino physics 
has undergone a revolution over the last decade or so by establishing the phenomenon of neutrino 
oscillations. The parameters relevant for three-flavour neutrino oscillations are three mixing an- 
gles, 612,613,623, and one CP-phase in the unitary lepton mixing matrix, and two mass-squared 
differences Am^i =m\ — m[ and Amgj =m^ — m\. In this note we summarize the determination 
of these parameters from present world neutrino oscillation data (Sec. ^, and we give a critical 
discussion of recent claims for possible hints in favour of a non-zero value of the mixing angle 613 
(Sec. |). 

2. Summary of three-flavour osciUation parameters 



Table I summarizes the results of two recent global fits to world neutrino data from Refs. [1,2]. 
For another recent analysis see [3]. 





Ref. [1] 


Ref. [2] (MINOS updated) 


parameter 


best fitibla 


3 a interval 


best fitibla 


3 a interval 


Am^i [lO^^eV^] 




7.05-8.34 


'•"'-0.21 


7.07-8.34 


Am]^ [lO^^eV^] 


±2.40^°;}? 


±(2.07-2.75) 


-2.39±0.12 
+2.49 ±0.12 


-(2.02-2.79) 
+(2.13-2.88) 


sin^ 012 


304+°"^^ 


0.25-0.37 


321+° °23 


0.26-0.40 


sin^ 023 


50+° °^ 

"•-'"-0.06 


0.36-0.67 


"•^'-0.06 


0.33-0.64 


sin^ 013 


QQ. +0.016 
"•"^-0.011 


< 0.056 


0.003 ±0.015 


< 0.049 



Table 1: Determination of three-flavour neutrino oscillation parameters from 2008 global data [1,2]. 



The latest data release from the KamLAND reactor experiment [4] has increased the exposure 
almost fourfold over previous results [5] to a total exposure of 2881 ton-yr and also systematic 
uncertainties have been improved. The Sudbury Neutrino Observatory (SNO) has released the 
data of its last phase, where the neutrons produced in the neutrino neutral current (NC) interaction 
with deuterium are detected mainly by an array of ^He NC detectors (NCD) [6]. The main impact 
of the new SNO data is due to the lower value for the observed CC/NC ratio, ((/>cc/0nc)^™ = 
0.301+0.033 [6], compared to the previous value (^cc/^c)'''" = 0.34 ±0.038 [7]. Since for 
neutrinos 0cc/^c ~ Pee ~ sin^ 0i2, adding the new data point on this ratio with the lower value 
leads to a stronger upper bound on sin^ 0i2. This explains the slightly lower best fit point for sin^ 0i2 
found in [1] compared to [2], since the latter does not (yet) include the SNO-NCD result. Data 
from SNO are combined with the global data on solar neutrinos [8] and with the recent results from 
Borexino [9], reporting a survival probability of the 0.862 MeV ^Be neutrinos of Pef'^'"^^ = 0.56 ± 
0.1. Fig. |I|(left) illustrates how the determination of the leading solar oscillation parameters 0i2 
and Am2i emerges from the complementarity of solar and reactor neutrinos. Spectral information 
from KamLAND data leads to an accurate determination of Atn^i with the remarkable precision of 
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Figure 1: Determination of the leading "solar" and "atmospheric" oscillation parameters [1]. We show 
allowed regions at 90% and 99.73% CL (2 dof) for solar and KamLAND (left), and atmospheric and MINOS 
(right), as well as the 99.73% CL regions for the respective combined analyses. 



8% at 3a, defined as (;c"PP'=' -;c'°*'='')/(-^"'''""' + -^ )■ KamLAND data start also to contribute to 
the lower bound on sin dn, whereas the upper bound is dominated by solar data, most importantly 
by the CC/NC solar neutrino rates measured by SNO. 

The MINOS experiment has reported new results on disappearance with a baseline of 
735 km based on a two-year exposure from the Fermilab NuMI beam corresponding to a total of 
3.36 X 10^° protons on target [10]. The data confirm the energy dependent disappearance of v^, 
showing significantly less events than expected in the case of no oscillations in the energy range 
< 6 GeV, whereas the high energy part of the spectrum is consistent with the no oscillation expec- 
tation. We combine the long -baseline accelerator data from MINOS as well as from K2K [11] with 
atmospheric neutrino measurements from Super- Kamiokande [12], using the results of Ref . [13]. 
In this analysis sub-leading effects of Am2i in atmospheric data are neglected, but effects of 613 
are included. Fig. [I|(right) illustrates how the determination of the leading atmospheric oscillation 
parameters 623 and lArngJ emerges from the complementarity of atmospheric and accelerator neu- 
trino data. The determination of | Anijj | is dominated by spectral data from the MINOS experiment, 
where the sign of \m\^ (i.e., the neutrino mass hierarchy) is undetermined by present data. The 
measurement of the mixing angle 623 is still largely dominated by atmospheric neutrino data from 
Super-Kamiokande with a best fit point close to maximal mixing. Small deviations from maximal 
mixing due to sub-leading three-flavour effects (not included in the analysis of [1, 13]) have been 
found in Refs. [2, 14], c.f. Tab. |I[ See, however, also Ref. [15] for a preliminary analysis of Super- 
Kamiokande. A comparison of these subtle effects can be found in Ref. [16]. At present deviations 
from maximality are not statistically significant. 

The third mixing angle 613 would characterize the magnitude of CP violation in neutrino os- 
cillations and is crucial also for the determination of the neutrino mass ordering. It is the main 
objective of upcoming reactor and accelerator experiments to measure this parameter. Fig. ^left) 
summarizes the information on du from present data, which emerges from an interplay of different 



3 



Three-flavour neutrino oscillation update 



Thomas Schwetz 



^ 4 

M 
> 

0) 



3 - 



<l 

1 



"T 1 — I — I I I I I I 1 1 — r 

30% 



90% CL (2 dof) 



global 



SK+K2K+MIN0S 



solar+KamL 
"+CHOOZ 

A. 



ZOr 



J I I I I II 1 1 1 



10' 



10 



Sin 



13 



0.15 



CD 

CM 



CO 



0.05 



1 1 1 
— f 


1 1 1 1 1 1 1 1 1 1 1 1 

— KamLAND 

\ global 






solar - 


1 1 1 


1., 1 


\ 

\ 

1 1 1 "^1^- 



0.1 



0.2 



0.3 

Sin e 



0.4 



0.5 



12 



Figure 2: Left: Constraints on sin 613 from the interplay of different parts of the global data. Right: 
Allowed regions in the (612 - 613) plane at 90% and 99.73% CL (2 dof) for solar and KamLAND, as well 
as the 99.73% CL region for the combined analysis. Am^i is fixed at its best fit point. The dot, star, and 
diamond indicate the best fit points of solar, KamLAND, and combined data, respectively. 



data sets. An important contribution to the bound comes, of course, from the CHOOZ reactor exper- 
iment [17] combined with the determination of lAmjJ from atmospheric and long-baseline exper- 
iments. The results on 613 [1] reported in Tab. |T] imply an upper bound of sin^ 613 < 0.035(0.056) 
at 90% CL (3a) and the best fit point is consistent with zero within 0.9a, while Ref. [2] finds the 
slightly stronger bound sin^ 613 < 0.049 at 3a and a Ax^ = 0.04 for 613 = 0. These differences can 
be understood by noting that SNO-NCD data is not included in the results from [2], see below for 
an explanation. 



3. Comments on possible hints for a non-zero value of 613 

A recent global analysis of neutrino data [18] obtains a hint for a non-zero value of dis at 
1.6a, which emerges from a 1.2a hint from solar-i-KamLAND data combined with a 0.9a hint 
from atmospheric-i-long-baseline-i-CHOOZ data. The hint from solar-i-KamLAND data can be un- 
derstood by the slight downward shift of the SNO CC/NC ratio due to the SNO-NCD data men- 
tioned in Sec. ^. From the combination of solar and KamLAND data we find a best fit value of 
sin^ 0,3 = 0.03 with Ax^ = 2.2 for 613 = which corresponds to a L5a effect (86% CL). We illus- 
trate the interplay of solar and KamLAND data in Fig. fright). The relevant survival probabilities 
are given by 

^ J cos^ 013 (1 — sin^20i2( sin^ 0)) solar, low energies / KamLAND 
cos 013 sin- 012 solar, high energies 



where = Am\^L/AE and (sin^ 0) w 1/2 for solar neutrinos. Eq. (3.1) implies an anti-correlation 
of sin^ 013 and sin^ 0i2 for KamLAND and low energy solar neutrinos. In contrast, for the high 
energy part of the spectrum, which undergoes the adiabatic MSW conversion inside the sun and 
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which is subject to the SNO CC/NC measurement, a positive correlation of sin^ and sin^ On 
emerges. As visible from Fig. ^(right) and as discussed already in [13, 19], this complementarity 
leads to a non-trivial constraint on 613 and it allows to understand the hint for a non-zero value 
of 013, which helps to reconcile the slightly different best fit points for 612 (as well as for Am^^, 
see Fig. |l]) for solar and KamLAND separately [18, 19, 20, 21]. This trend was visible already in 
pre-SNO-NCD data, though with a significance of only 0.8a, see v6 of [13]. 

The hint arising from atmospheric data is, in our opinion, much more controversial. This 
preference for a non-zero 613 value was first noted in Ref. [14], and is illustrated in Fig. 24 of 
that paper (see also [22]). Our attempts to reproduce that figure with our own numerical codes are 
shown in the left panels of Fig. |3[ We have performed three different kinds of calculations: 

(A) old rates, old (black line). This case follows closely the analysis presented in Ref. [13], 
except for the fact that here we have also included sub-leading Am\^ effects. Note that the 
data samples and the definition used here are the same as in Ref. [14]; in particular, SK-I 
data are divided into 55 energy and zenith bins. 

(B) new rates, old (red line). This case differs from the previous one for the details of the 
event rate calculations, which are now performed as described in the appendix of Ref. [2] ; 
on the other hand, the x^ definition and the number of data points are unchanged. 

(C) new rates, new x^ (blue line). This case fully coincides with the three-neutrino analysis 
presented in Ref. [2]. The treatment of the theoretical uncertainties in the x^ definition differs 
considerably from the previous two cases, and the number of data points is now increased to 
90. 

For definiteness we consider only normal hierarchy, and we have explicitly verified that the in- 
clusion of K2K and/or MINOS does not affect the results, compare also [14] (without MINOS) 
and [18] (with MINOS). Our analysis shows that: 

• the results of case (A) are qualitatively very similar to those shown in Fig. 24 of Ref. [14]. 
In particular, the data favor cos 5 = — 1 and disfavor cos 5 = +1, indicating that a relevant 
contribution to the effect comes from the interference term A3 (see Eq. 44 of Ref. [14]). 
Indeed, from the lower-left panel we see that the hint for non-zero 613 disappears for Km\^ = 
0, in agreement with [14]; 

• on the other hand, the strength of the preference for non-zero 613 in our analysis (A) is 
only 0.5a, much weaker than the 0.9a found in Ref. [14]. Since the definition and the 
experimental data are the same for both analyses, we conclude that the discrepancy must 
arise from differences in the rate calculations. Indeed, a comparison between our cases (A) 
and (B) in Fig. ^ shows that the statistical relevance of the signal strongly depends on the 
details of the Monte-Carlo: for example, for case (B) the hint of non-zero 613 is only at the 
0.2a level, i.e. practically inexistent; 

• details of the data binning and the theoretical uncertainties are also very important, as can be 
seen by comparing cases (B) and (C), which are based on identical rate calculations. In case 
(C) the preference for non-zero 613 is again as strong as in case (A), but now the hint persists 
also for A/M21 = 0. 
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Figure 3: Zenith distribution for muhi-GeV e-like events (upper panels), and Ax^ dependence on 
cos 5 sin 013 for Amji = 8 x 10^^ eV^ (middle panels) and for Amjj = (lower panels). The fits include 
the CHOOZ data as well as the full atmospheric data samples for SK-I (left panels), SK-I with modified 
multi-GeV e-like data (central panels) and SK-I H- SK-II (right panels). The black, red and blue lines corre- 
spond to different analyses. See text for details. 

It is quite difficult to trace the preference for non-zero d\3 to the effect of a single data sample. 
Clearly, the asymmetry between the two values cos 5 = ±1 shows that the interference term A3, 
mostly relevant for sub-GeV data, plays an important role. On the other hand, other samples are 
affected by 613, and as we have seen the details of the statistical analysis determine whether the hint 
appears also when the interference term is suppressed {e.g., for Am^i = 0), or not. As an example, 
in the upper panels of Fig. ^ we show the contribution of multi-GeV e-like data to the effect. It is 
clear from this plot that the first two bins of this sample are better fitted with a non-zero value of 
013. To assess the relevance of this contribution, in the middle and right panels of Fig. || we repeat 
the calculations described above for a different choice of experimental data: 
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• in the middle panels, we artificially reduce the first two bins of SK-I muIti-GeV e-like events 
by 10%. As can be seen, with this simple change the hint in favor of non-zero completely 
disappears, irrespective of the details of the rate calculations and of the statistical analysis, 
and also of the inclusion of sub-leading Am^i effects. We conclude therefore that the excess 
in these two bins somewhat "triggers" the hint for non-zero 613; 

• in the right panels, we repeat our fits for the combination of SK-I and SK-II data. As can be 
seen from the upper panel, the "excess" in the first bins of multi-GeV e-like data no longer 
appears once SK-II data are also included. The lower panels show that the results are very 
similar to the "fake" SK-I data displayed in the middle panels, and in particular no preference 
for non-zero 613 is present. 

In the light of these results, we conclude that (1) the statistical relevance of the claimed hint for 
non-zero 613 from atmospheric data depends strongly on the details of the rate calculations and of 
the analysis, and (2) the hint is the result of a statistical fluctuation in SK-I data which disappears 
after the inclusion of SK-II data. Since in the atmospheric neutrino analysis used in [1, 13] no hint 
for 613 > is present, the indication from solar data becomes weakened in the global analysis, 
leading to the du result quoted in Tab. |I|, consistent with zero at 0.9a. We beheve that the present 
"indication" for a non-zero 613 should not be taken more serious than a ~ la effect (which of 
course implies that it will be true with a probability of 68%) and has to wait for being confirmed 
or refuted by upcoming experiments. 
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